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ACOUSTIC  INSTABILITIES  AS  A  SOURCE  OF  TURBULENCE 


IN  REDUCED  DENSITY  CHANNELS 

I.  Introduction 

In  experiments  at  NRL  hot,  reduced  density  channels  have  been  formed  in  air  and  argon  by 
ohmic  heating1  and  in  nitrogen  by  laser  heating.2  Nominally  such  channels  could  persist,  in  pressure 
equilibrium  with  the  cold  surrounding  gas,  until  they  cooled  by  heat  conduction.  In  fact  they  become 
turbulent  and  dissipate  —2000  times  faster  than  heat  conduction  would  dictate.  Boris  and  Picone3  have 
shown  that  spatial  asymmetries  in  the  initial  heating  of  the  channels  can  cause  such  turbulence.  How¬ 
ever  the  appearance  of  turbulence  in  channels  with  relatively  symmetric  initial  healing  suggests  the  pos¬ 
sibility  of  additional  sources. 

Electrical  measurements  on  the  channels  in  air  show  that  they  do  not  decay  to  an  equilibrium 
electron  density  for  —100  ^ s .  In  argon  significant  amounts  of  energy  can  be  trapped  (resonant  self 
trapping)  and  concentrated  in  metastable  levels.  Air  and  nitrogen  can  also  store  energy  in  vibrational 
states  with  long  relaxation  times.  It  has  been  recognized  that  acoustic  wave  growth  can  occur  in  a  gas 
maintained  in  a  nonequilibrium  state.4  This  paper  demonstrates  that  such  growth  can  occur  in  a  gas 
relaxing  from  a  nonequilibrium  condition.  The  theory  is  then  applied  to  the  channel  situation  to  see  if 
energy  "frozen  in"  during  the  channel  expansion  is  sufficient  to  cause  a  disturbance  to  amplify  to  the 
point  of  becoming  turbulent. 

Finally  we  show  the  channels  are  stable,  if  the  gas  remains  in  thermal  equilibrium,  despite  the  gas 

possessing  a  temperature  dependent  specific  heat. 

Manuscript  submitted  February  5,  1982. 
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II.  A  Description  of  the  Acoustic  Instability 

As  a  simple  model  which  will  demonstrate  the  growth  phenomena  we  consider  a  diatomic  gas  pos¬ 
sessing  a  long  lived  vibrational  level.  We  suppose  this  level  may  be  depopulated  only  by  binary  colli¬ 
sions  whereupon  its  energy  is  given  up  to  the  internal  kinetic  energy.  If  the  gas  properties  remain  spa¬ 
tially  uniform  a  slow  relaxation  to  thermal  equilibrium  occurs.  A  sound  wave  passing  through  the  gas 
causes  local  perturbations  in  the  temperature  and  density  and  thus  local  perturbations  in  the  relaxation 
rate.  A  region  of  compression  speeds  the  relaxation  rate  delivering  an  extra  energy  kick  to  the  subse¬ 
quent  re-expansion.  On  the  other  hand  slowing  of  the  relaxation  rate  in  the  rarefactions  assists  their 
re-collapse.  The  sound  wave  thus  grows  by  diverting  some  of  the  energy  from  the  decaying  vibrational 
state  into  bulk  motion  and  away  from  internal  kinetic  energy. 

Sound  waves  are  normally  damped  by  viscosity  and  by  heat  conduction  between  the  hot, 
compressed  and  cold,  rarefied  zones.  We  will  neglect  viscosity  and  heat  conduction  in  the  present 
analysis.  This  approximation  must  thus  be  justified,  in  any  application,  by  comparison  of  the  predicted 
growth  rate  to  the  so-called  classical  absorption.' 

Finite  relaxation  times  for  internal  degrees  of  freedom  normally  introduce  additional  damping  due 
to  the  loss  and  retrieval  of  energy  to  and  from  these  states.  This  effect  is  outweighed  in  the  present 
case  where  the  vibrational  level  is  too  high  lying  to  be  accessible  at  the  post  expansion  channel  tem¬ 
perature. 

III.  Linear  Stability  Analysis 

We  describe  the  behavior  of  the  gas  by  the  usual  Eulerian  variables  density,  p,  pressure,  p,  and 
velocity,  v,  plus  a  specific  internal  energy  exclusive  of  the  vibrational  level,  «.  and  a  specific  energy  in 
the  vibrational  level  alone.  «  *  The  latter  energy  is  related  to  the  density  of  excited  atoms,  n *.  via 

t i"!  (1) 

P 

where  £is  the  excitation  energy.  The  rate  equation  determining  n‘ 
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Bn* 

Bt 

may  be  expressed  in  terms  of  «  *,  p,  and  «.  (We  neglect  the  last  term  representing  relaxation  damping 


— (T  n'nv  +  it  *n(n  —  n*)v 


(2) 


since  we  assume  <r  *  «  <r  and  n*  >  n.) 


Be*  ,  /  32A0t 

Bt  <T€  P  "v  5irA  m 

where 

<t  *  =*  cross  section  for  collisional  excitation 
<t  =**  cross  section  for  collisional  de-excitation 
n  ==  total  number  density 
v  =»  mean  velocity  between  identical 

molecules  = 

A0  »  Avogadro’s  number 
A  »  gram  molecular  weight 
m  =»  mass  per  molecule 

Equation  (3)  may  be  rewritten  for  a  moving  fluid  by  following  standard  procedures6 


De* 

Dt 


A  e  *  pV7 


where  D/ Dt  is  the  material  derivative  given  by 


(rate  equation) 


and  we  have  defined 
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Other  governing  equations  in  addition  to  (4)  include 


and 


-  -p  V  v  (conservation  of  mass) 
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Dv 

D/ 


-V  p  (conservation  of  momentum) 
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where  viscosity  has  been  neglected.  The  conservation  of  energy  equation  contains  an  extra  term  to 
retleci  the  transfer  of  energy  from  the  vibrational  level  into  internal  kinetic  energy 

p  —  — />V  •  v  -  (>  ~p~  (conservation  of  energy)  (7) 

and  heat  conduction  has  been  neglected. 


One  additional  equation  is  necessary  to  provide  a  closed  set.  Since  the  vibrational  energy  interacts 
with  the  fluid  only  via  the  decay  represented  by  the  last  term  in  (7),  the  fluid  variables  are  interrelated 
as  in  an  equilibrium  gas.  We  may  therefore  adopt  the  usual  equation  of  state 


e 


(equation  of  state). 


(8) 


A  slowly  relaxing,  spatially  uniform  medium  represents  a  succession  of  equilibrium  states  to  a 
sound  wave  whose  period  is  short  compared  to  the  decay  time.  The  fluid  variables  may  be  represented 
as  a  slowly  changing  quasi-equilibrium  value  plus  a  rapidly  changing  perturbation 


«*  =  «o*+«r 

etc.  (9) 

Two  possible  types  of  quasi-equilibria  are  of  interest,  constant  pressure,  />,>.  and  constant  density.  p„. 

The  former  more  nearly  represents  the  situation  within  the  channel  but  either  is  acceptable  under  the 

slow  relaxation  approximation.  The  latter  leads  to  simpler  algebra  and  we  shall  develop  that  case.  Our 

assumption  of  a  constant  density,  spatially  uniform  equilibrium  implies 


$1-0  <10, 

dr 

Vp„-  Vt„=  V«„*-  0.  (11) 

These  assumptions,  applied  to  the  governing  equations.  (4)  through  (8).  show  that  the  equilibrium  is 


further  described  by 
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(15) 

(16) 


If  we  substitute  the  split  forms  of  the  variables,  (9),  into  the  governing  equations,  (4)  through 
(8).  eliminate  those  terms  which  describe  the  equilibrium  behavior  by  invoking  (10)  through  (16),  and 
retain  only  terms  to  first  order  in  the  perturbations,  we  are  lead  to  the  linearized  equations  for  the  per¬ 
turbation  variables 


foi' 
8 1 


=  -  A 


ml -  ,  •  / -  ,  g  OP  0  €  1 

Pif- 1  v€o  +  foPiv€o  +  2yfe~ 


^Pt  _ 

8  v , 

p017  =  -Vp, 

8€  |  8e  * 

Po^T=-PoV  — 


(17) 

(18) 

(19) 

(20) 
(21) 


Pj_  _  P_i_ 

Po  P  o 

We  need  only  consider  planar  disturbances  to  demonstrate  the  growth  phenomena.  Given  linear  equa¬ 
tions  it  is  sufficient  to  consider  separately  each  Fourier  component  of  any  sonic  disturbance.  Substitut¬ 
ing  a  normal  mode  form  for  each  perturbation  variable 


into  (17)  through  (21)  yields 


«,***  e  ,*((«. k)e'ikx-'“'' 
etc. 
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(22) 
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(24) 

(25) 

(26) 

(27) 


The  perturbation  variables  are  now  understood  to  be  functions  of  <o  and  k. 
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Equations  (23)  through  (27)  may  be  solved  for  any  perturbation  variable,  the  pressure  proves  to 
be  a  convenient  choice.  We  solve  (23)  for  e  *,  and  use  that  result,  and  (27),  to  eliminate  « (  and  «  *  in 
(26).  Equations  (24)  and  (25)  may  be  used  to  eliminate  pi  and  v,  resulting  in 

UJ'Poeo  ,  6-4eoV*o  At  i  ,  Po-4  e»k~yje Tt 

Pi  -  +  - - ; - - e 0Ar- - Ar  4-  — - j=-  =  0  .  (28) 

Po  2^q(/oi  •‘ipo’v^o)  Po  2 — 

Setting  the  bracketed  term  to  zero  yields  the  dispersion  relation  between  to  and  k.  This  becomes,  upon 
further  simplification 

<o2(2/toV«o  -  2« ,).4 p i)  +p0.4«l,*)  -  fc:(2/c,Ja»v'«t)  “  2c~Ap,£ c|  -  4/>(<e„*)  -  0  (29) 

where  we  have  introduced  the  equilibrium  sound  speed 

'  Po  Pit]  y  Pn  ,, 

t  —  -  «o  +  —  =  -  (30) 

Po*n  Pu  Pi 

(Notice  that  turning  off  the  coupling  between  «  and  « *,  by  setting  A  to  zero  recovers  the  dispersion 
relation  for  normal  sound  waves,  to2  =  c:k'.) 


The  solution  of  equation  (29)  requires  to  to  be  complex.  We  may  substitute  to  =  to,  +  it o,  in  (29)  and 
separately  equate  the  real  and  imaginary  parts  of  the  resulting  expression  to  zero.  This  gives 

to  •  ( 6<o  ,-JT~q  +■  2e  o  .4  p  o  —  P  1)  4*0*) 

-  A 2 1 2  r2<o . v'«Ti  +  2c:  4p.fi,)  +  .4p,ie,i*) 

—  to  - ( 2to ,-y/e ,j  +  2  e,j-4p o  —  Po  4  «,)")  =0  (31 ) 


3<o  ,2v«o  +  ( 2e  o  4  p  o  -  p,>.4€n*)co,  +  yje„(c:k2  -  to2)  =  0. 
Equation  (31)  may  be  used  to  eliminate  to2  in  (32)  with  the  result  that 


where 


<0 4-  V/ oj  ~  4-  .Vcu,  4-  Q  =  0 


M  =■  ^po-v/*rt  ~ 


'4po*o 


V  =  4~  _  P.i  4-tn  +  Pn  4~  en  ‘  C2A2 

4  4  +  I6«n  +  4 


<?=  - 
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The  constant  term,  Q ,  will  always  be  negative,  thus  there  will  always  be  a  positive  value  of  <a, 
corresponding  to  wave  growth. 

IV.  Application  to  Laser  Formed  Channels  in  Nitrogen 


'The  de-excitation  cross  sect  ion.  rr  .  introduced  in  Section  III  and  p  m  arc  related  in  via  <r  =  Pm IT  mil- 
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The  channel  expansion  takes  —  30  m  sec  so  vibrational  levels  may  first  become  excited  and  then  frozen. 


In  nitrogen  only  the  lowest  vibrational  level  (E/k  —  3340 A  )  will  be  appreciably  excited  under  the 
initial  conditions.  As  an  upper  limit  on  the  possible  vibrational  energy  we  calculate  the  equilibrium 
population  for  this  state 


If  this  fraction  of  the  molecules  remains  excited  during  the  expansion  we  find 


•  ",  ,  -r  ■'() 

fo  =*  —  kT,  — 
n,  A 


1.5  x  10* 

g 


W'e  now  calculate  the  growth  rate.  Under  these  conditions  only  the  linear  term  in  <331  and 
the  last  term  in  (35)  are  significant.  The  growth  rate  becomes  independent  of  wave  number,  k,  and  is 
given  by 


where  we  have  used  (see  Eq.  12) 


to  evaluate  I  . 


I  Poly  -  n  «i* 

- — - =  36  sec 

4  VV« o 


r  P  o  V  6  a 


This  growth  rate  is  clearly  too  small  to  have  any  effect  over  the  300  ms  during  which  the  nitrogen 


V.  Stability  of  Gases  with  Temperature  Dependent  Specific  Heats 


For  low  frequencies  or  fast  relaxation  rates  the  gas  remains  in  thermal  equilibrium  as  it  adjusts  to 
the  disturbance.  We  may  include  the  contributions  from  all  the  internal  degrees  of  freedom  in  e  and 
set  e  *  to  zero.  The  specific  energy  now  acquires  a  non  linear  temperature  dependency 


« oc  r 
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where  a  >  1.  If  we  combine  (43)  with  the  gas  law,  p  =  p  RT,  and  linearize  the  resulting  expression  we 
find 

Pi_  _  Pj_ 

Po  Pa 

The  dispersion  relation  for  perturbations  from  fixed  equilibrium  is  then  given  by  substituting 
«o“  (lf  o  and«,*=  0  in  (29).  Ordinary  sound  waves  now  result  (tu  =  c/c)  with  a  modified  sound  speed 


(44) 


j  Pa 

J  Pa 

J  a  Pit  a 

fje„  a  — 
Pa 

The  medium  is  now  stable  to  any  planar  disturbance. 


(45) 


The  energy  principle1*  permits  us  to  demonstrate  stability  in  complex  geometric  situations  without 
having  to  solve  the  equations  of  motion.  We  shall  apply  it  to  a  column  of  hot  gas.  with  a  temperature 
dependent  specific  heat,  embedded  in  a  cold  surrounding  medium. 


The  Lagrangian  density  in  a  region  of  uniform  pressure  and  density  may  be  expressed  as10 

* 

L  -  T-  V-  \Pnv2  +  |v0  .  (V  -t,)2  (46) 

2  2  av  0 

where  y  is  a  small  fluid  displacement  from  equilibrium  and  i'()  =  l/p„  is  the  specific  volume.  Proceed¬ 
ing  as  in  the  derivation  of  (44)  we  find 

2  dp,  dv 

=  a«o  — ~  4 -  (47) 

P  o  ro 

But  having  neglected  heat  conduction,  the  first  law  of  thermodynamics  states 


Equating  (47)  and  (48)  gives 


so  that 


d«  =  -  Podv. 


ftp  Pa  .  a  eo 

- - Pa  + - 

bv  a  f  o  v0 


(48) 


(49) 
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1  ,2  1  P0V0 


r  1  ‘l  1 

L  -  -r  Po»I  “  T 

2  2  a  c0 


Po  + 


a  €0 


(V  •  V2. 


(50) 


Applying  Lagrange's  equation  for  continuous  media"  to  (50)  we  recover  the  wave  equation  with 
the  sound  speed  given  by  (45). 

The  energy  principle  states  that  the  motion  of  the  system  will  be  bounded  if  the  total  potential 
energy  for  the  system  increases  for  any  disturbance.  The  potential  energy  is 


t  olumn  h 


Po  + 


(V  •  T,)2 
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i///A 


Po1', 

a 


Po  + 


a  e,. 


(V  ■  Tf)2 


(51) 


where  the  subscripts  /;  and  c  refer  to  hot  and  cold  respectively.  The  potential  energy  is  positive  for  any 
compressive  (i.e.  sonic)  disturbance  and  zero  for  any  incompressible  flow.  These  situations  correspond 
to  bounded  oscillations  and  persistent  but  non  amplifying  flows. 


VI.  Conclusions 

An  excess  of  nonequilibrium  internal  energy  in  a  gas  transforms  the  relaxational  damping  of 
sound  waves  into  a  growth  phenomena.  A  model  diatomic  gas  has  been  considered  which  possesses 
one  vibrational  level.  The  small  signal  growth  rate  for  normal  modes  (plane  waves)  has  been  calculated 
from  the  rate  equation  and  the  equations  of  motion.  The  growth  proves  to  be  slight  for  the  conditions 
that  obtain  within  the  laser  heated  reduced  density  channels.  Thus  it  appears  that  air  chemistry  plays 
no  role  in  creating  the  turbulence  observed  in  these  channels.  Meiasiable  energies  an  order  ol  magni¬ 
tude  or  more  larger  than  the  kinetic  energy  are  required  to  produce  a  significant  efTect. 

An  energy  principle  analysis  shows  that  the  internal  energy  stored  in  rapidly  relaxing  states  cannot 
cause  instabilities  i.e.  a  gas  which  remains  in  thermal  equilibrium  but  possesses  a  temperature 
dependent  specific  heat  does  not  constitute  an  amplifying  medium  for  sound  waves. 
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